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The properties of nuclei embedded in an electron gas are studied within the relativistic mean- 
field approach. These studies are relevant for nuclear properties in astrophysical environments 
such as neutron-star crusts and supernova explosions. The electron gas is treated as a constant 
background in the Wigner-Seitz cell approximation. We investigate the stability of nuclei with 
respect to a and (3 decay. Furthermore, the influence of the electronic background on spontaneous 
fission of heavy and superheavy nuclei is analyzed. We find that the presence of the electrons 
leads to stabilizing effects for both a decay and spontaneous fission for high electron densities. 
Furthermore, the screening effect shifts the proton dripline to more proton-rich nuclei, and the 
stability line with respect to /3-decay is shifted to more neutron-rich nuclei. Implications for the 
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creation and survival of very heavy nuclear systems are discussed. 
PACS numbers: 21.10.Dr, 21.10.Ft, 21.10.Tg, 23.60.+e, 26.50,+x, 



*Elcctronic address: buervenich@fias.uni-frankfurt.de 
^Electronic address: mishustin@fias.uni-frankfurt.de 
■^Electronic address: greiner@fias.uni-frankfurt.de 



1 



I. INTRODUCTION 



Nuclear astrophysics lies at the heart of understanding the origin of the physical world 
and our mere existence. One of its key questions is 'How and where did the chemical 
elements originate?'. Many challenging questions concerning the processes and mechanisms 
for nucleosynthesis are in the focus of modern research. In this respect, supernovae and 
neutron stars are the most interesting objects to study. 

The first attempt to explain abundances of chemical elements in the solar system was 
made by Geoge Gamow in 1946-48 (see e.g. Ref. jj). He proposed a mechanism that was 
based on a continuous building-up of chemical elements by neutron capture in the early 
universe. Later it became clear that this mechanism cannot explain the production of heavy 
elements because of very high specific entropy of the early universe, S/B~ 10 10 , where S is 
the total entropy and B is the net baryon number. The present understanding of the element 
creation says that the lightest elements (up to He, and, partly, Li) were formed during the 
first moments of the universe expansion, immediately after the Big Bang. Heavier elements 
up to oxygen are predominantly produced in thermonuclear reactions in stars like our Sun. 
Elements up to iron and nickel could be produced in more massive stars. The origin of 
heavier elements such as gold and uranium remains in the focus of current research. It is 
widely believed that heavy elements were mostly synthesized in the course of supernova 
explosions in the so-called r-process (subsequent neutron capture by stable and unstable 
nuclei, see a recent review Q). 

A type II supernova explosion is one of the most spectacular events in astrophysics, with 
huge energy release of about 10 53 erg or several tens of MeV per nucleon |3j. When the core 
of a massive star collapses, it reaches densities several times larger than the normal nuclear 
density po — 0.16 fm~ 3 . The repulsive nucleon-nucleon interaction gives rise to a bounce-off 
of matter and formation of a shock wave propagating through the in-falling stellar material, 
predominantly Fe. 

Hydrodynamical simulations (see e.g. Refs. {4], 5]) show that during the collapse and 
subsequent explosion the temperatures T ps (0.5 4- 10) MeV and baryon densities ps ~ 
(10~ 5 4 2)p can be reached. Unfortunately, these simulations do not produce successful 
explosions yet, even when neutrino heating and convection effects are included. This means 
that some important physics of this process is still missing. 
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Besides supernova explosions, proto-neutron and neutron stars are interesting objects 
from a nuclear physics point of view. They are quite similar to huge, extremely neutron-rich 
macroscopic nuclei that gain additional stability and binding due to the gravitational force. 
Neutron stars are very compact objects with a central density of about p = 10 15 g/cm 3 , 
a typical radius of R = 10 km and masses up to 2 solar-masses. Proto-neutron stars are 
newly-born neutron stars formed in the course of a supernova explosion. They are somewhat 
bigger than neutron stars and have temperatures up to 30 MeV. The regime of interest for 
our considerations are baryon densities in the range 0.001 — 0.5po, where very heavy and 
neutron-rich nuclei may be present {(], Q, [§], [9]. Nuclear pasta phases, i.e., nuclear matter 



in various geometries such as slabs and parallel plates, reminescent of pasta [lOl . I 111 . Il2l . 1 1 31 ] 
start at densites slightly above 0.5 po- I n this regime the neutrons need to be taken into 
account. 

Since the heaviest elements occuring in nature are Uranium isotopes, there should have 
been corresponding extreme conditions for their creation and persistence, in supernova ex- 
plosions, crusts of neutron stars, or other sites. Thus we may wonder if very heavy and 
superheavy systems can and will be produced under such conditions too. How are the 
properties of such nuclear systems altered in dense environments? 

The properties of nuclei in astrophysical environments have enjoyed continuous interest 
for over more than 30 years. The seminal work by Negele and Vautherin [ijl layed out the 
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or microscopic studies of nuclei in stellar environments, see e.g. Refs. 
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Q, 0]- The investigation of the rich properties of nuclear systems under extreme 



conditions in astrophysical environments is on the rising path. Since many astrophysical 
calculations, from r-process to dynamical simulations of supernova explosions, depend on 
theoretical nuclear input, it is important to understand how nuclei subject to these special 
conditions differ from nuclei studied in the laboratory on earth. 

In this paper we focus on the physics of nuclei embedded in a dense electron gas. This 
investigation is important by itself as well as for further studies including the neutron gas. 
We consider nuclei across the perodic chart up to superheavy nuclei embedded in a Wigner- 
Seitz cell with constant electron density. The presence of the electrons will effect the location 
of the /5-stability line and the proton dripline as well as decay modes such as a-decay and 
spontaneous fission. 

The paper is structured as follows: In Section [Til we present the relativistic mean-field 
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model which is employed in the calculations as well as the concept of the Wigner-Seitz cell 
and its concrete implementation for spherical and deformed nuclear systems. In Section IHIt 
we discuss the /3-equilibrium condition and present results for the /3-stability line and the 
proton and neutron drip lines in the presence of electrons. The evolution of the a-decay mode 
in heavy nuclei as a function of the electron Fermi momentum is demonstrated in Section 
IIV Al In Section HVBt deformed nuclei are considered and spontaneous fission under the 
influcence of the electron background is studied. Finally, in Section |V] we conclude and give 
an outline of future research directions. 



II. THE FRAMEWORK 



A. The RMF model 



Over the years self-consistent mean-field models have reached high predictive power. They 
can be applied from medium-light systems up to superheavy nuclei and to systems ranging 
from the proton drip line to the neutron drip line. They are based on the formulation of an 
effective nucleon-nucleon interaction meant to be employed in the Hartree or Hartree-Fock 
treatment of nuclear systems. The modern way of formulating them is in terms of an energy 
functional that can incorporate terms which can not be constructed via a two- (or three- or 
four-) body force. 



In this paper we employ the relativistic mean-field (RMF) model [2ll . |22| |. The effective 
in-medium nucleon-nucleon interaction is parametrized via the exchange of several meson 
fields: scalar-isoscalar (cr), vector- isovector {uj^) and vector-isovector (p M ). 

This model is based on an effective Lagrangian of the form 

a 
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The field tensor for the u> meson is defined as = d^uj u — d u uj^, and similar definitions hold 
for the field tensors of the p meson and the photon. The parameters g a , g^, g p , m a , b, c are 
fitted to experimental data of nuclear ground-state observables. The masses of the p and to 
mesons are fixed at the experimental values, since the performance of the model is not quite 
sensitive to their values. As it is written down in the Hartree approximation, the meson 
fields are treated as classical potentials and the nucleons are represented by Dirac spinors. 
The isoscalar-scalar a meson delivers the intermediate-range attraction, while the isoscalar- 
vector uj meson is responsible for the short-range repulsion. The isovector-vector p meson 
couples to the isovector nucleon density and thus parametrizes the isovector properties of 
the model. The photon is coupled in standard fashion. For all boson fields, no (explicit) 
exchange terms are taken into account. 

It is worth mentioning that these meson fields have only loose correspondence with the 
physical meson spectrum. Mean-field models employing contact interactions between nu- 
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states 



rave a comparable predictive power for nuclear ground-state observables and excited 
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The w a denote occupation probabilities of the nucleon states and originate from the 
treatment of pairing. We employ BCS pairing with a density-independent 5-force, see Ref. 



261 ] for details. In time-reversal even-even systems, only the time components of the vector 
mesons are non-zero. Proton and neutron states are not allowed to mix, hence only the 
third components of the isovector-vector p field associated with t 3 survives. 
The single-particle equation for the nucleons reads 

[17 • d + m N + g a o + gjj ^ + g P p° 3 r 3 + eA° — ^70]^ = Hoip (2-2) 

The binding properties of nuclear matter and nuclei are generated from the strong scalar 
and vector fields, Us = g^ ~ —350 MeV, Uy = g u oj Q ~ +300 MeV, which add up to a 
normal nucleon potential Un = Us + Uy ~ —50 MeV. They add up with the same sign 
to generate the strong spin-orbit potential in nuclei, which (in the nonrelativistic limit) is 
given by 

V ls (x^-(Vs-V vl0 )t-s (2.3) 
ar 

This spin-orbit force emerges from the covariant formalism with the right sign and magni- 
tude without introducing additional parameters. This is an important consequence of the 
relativistic description. 



In order to account for the electron background, we modify the source term of the photon 
field by adding the electron density, i.e. 



= -e(p p - p e ) = -ep ch , (2.4) 

which is the Poisson equation for the electrostatic potential <fi = Aq. Further down we 
discuss the concrete implementation in different symmetries. All other equations remain the 
same as in vacuum. 

Calculations are performed in coordinate space, but the derivatives are performed as ma- 
trix multiplications in Fourier space. We employ the mean-field parametrization NL3 in these 



studies 27] which delivers accurate values for nuclear ground-state properties. Furthermore, 
since the external electron background couples only electromagnetically to nucleons, no ad- 
ditional parameters need to be introduced, and no readjustment of the present parameters 
is needed. 



B. The Wigner-Seitz approximation 

In this paper we study properties of nuclei embedded in a uniform electron gas at zero 
temperature. The calculations are performed within the Wigner-Seitz (WS) approximation 
by dividing the system into WS cells, each containing only one nucleus and the number of 
electrons equal to the nuclear charge Z, see Fig. [I] for an illustration. 

The construction of the WS cell is aimed at an approximate and convenient segmentation 
of an ensemble of nuclei surrounded by the uniform background of electrons. Its construction 
should ensure that the physics contained in the Wigner-Seitz cell is to an optimally large 
degree isolated from the surrounding. This requires, for example, charge neutrality of the 
cell. However, there could be also other requirements which minimize the electrostatic inter- 
actions with neighbouring cells (see below). Recently the validity of the WS approximation 



for the inner neutron-star crust has been discussed, e. g. in Ref. 

At electron densities considered in this work, the Coulomb potential due to the nucleus 
is a small perturbation compared to the Fermi energy of the electrons. Thus the electrons 
will only be weakly affected by the presence of the nucleus which justifies the approximation 
of constant density. 

In the numerical implementation, however, we use a Fermi-type distribution of the elec- 
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FIG. 1: (COLOR ONLINE): Illustration of a Wigner Seitz cell for a spherical nucleus (left) and 
an axially deformed nucleus (right). The inner dark filled region denotes the charge density of 
the nucleus, the outer region denotes the extension of the electron background surrounding the 
nucleus. The figure is not to scale. 

tron density with a smooth boundary. The reason is that on a spherical grid in coordinate 
space, the cell radius cannot be fixed exactly (it is limited to any value that is a multiple 
of the grid spacing Ar). A smooth surface region allows us to realize the electron shell for 
any grid spacing and radius. The parameters of the shell are chosen to satisfy the charge 
neutrality condition J dV p c h(r) = 0, where p c h(r) is the total charge density. Since in the 
computer code the derivatives are calculated in Fourier space, the numerical realization of 
the smooth surface region is not very sensitive to the chosen grid-spacing. This construction 
of the electron density allows us to calculate nuclei throughout the nuclear chart without 
introducing artifical jumps in the energy caused by a step-like change of the WS cell radius. 

1. The spherical cell 

For a given nuclear charge Z the WS cell radius Rc is uniquely determined by the charge 
neutrality requirement 



Expressing the electron density in terms of the electron Fermi momentum kp, p e = -^k^ 
one can write the cell radius as 




(2.5) 




(2.6) 
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The parametrization with a smooth surface, as used in this work, is given by 

Pe{r) = Pe \ r _ rc) , (2.7) 
1 + exp v - ; 

where p e Q is the background electron density, tq is adjusted in order to fullfill charge neu- 
trality, and a is the diffuseness parameter of the distribution, taken to be 0.45 fm. The 
parameter a is chosen such that the charge density is dropping from 90% to 10% of its value 
at r = within the diffuseness interval o = 4.4* a = 2 fm. We have checked that the results, 
in particular the differences of nuclear binding energies, are not sensitive to the choice of a. 
The Poisson equation (12. 4p has to be solved numerically. 

While these calculations can be performed quite easily for high electron densities, for low 
electron densities when the WS cell radius becomes quite large, numerical problems arise. 
For example, for a tin nucleus, at kp = 0.1 fm -1 , the WS cell radius is Rc ~ 70 fm, while 
for kp = 0.01 fm -1 it corresponds to Rc ~ 700 fm. 

An alternative way of calculating the effect of the electron gas is to add the potential of 
a homogenously charged sphere of electrons to the electric potential of the proton charge 
distribution, i.e., first calculate the proton electric potential, 

A0 P = -ep p (2.8) 

then add the potential due to electrons, 

(f> = <j>p + 0e (2.9) 

where cj) p is the solution of the Poisson equation for the proton charge density. The electric 
potential caused by the uniform background of electrons reads 

-Ze [3 1 r 



e (r) = { ^ UCLZ ZK c^ - (2.10) 

lf> r>Rc 

As can be seen from Fig. [2j the two different ways of calculating the total electric 
potential, Eq. H2.4[) and Eq. (12.91) . lead to identical electric potentials (apart from the 
correction due to the smooth surface region in Eq. (12.41) ) and, therefore, to the same 
structure of the nucleus. This demonstrates consistency and accuracy of the treatment of 
the electron background in our numerical calculations. 

It is clear that the total energy of the cell should include the interaction energy of the 
electrons with the electromagnetic field, which does not appear when the electrons are 
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FIG. 2: The electrostatic potential 4> shown for the doubly-magic lead nucleus for the cases of no 
electrons, for the inclusion of the electron density (p) and for the inclusion of the external electron 
potential {<f> e ), respectively, at kp = 0.1 fm -1 . 



employed as an external potential, Eq. (12.91) . This interaction energy can be easily calculated 
for the case when the (spherical) electron and proton charge distributions are replaced by 
step-functions: 

Pe(r) =p e0 O(Rc-r), p p (r) = p p0 6(R A - r) (2.11) 

where Rc and Ra are the cell radius and nuclear radius, respectively. An elementary calcu- 
lation yields for the total electrostatic energy of the cell 



E, 



3 Z 2 e 2 



Coul 



c(^), c(*) = l 



3 1/3 1 

-x 1 + -x 



(2.12) 



5 47i\R A >po' 

This expression has the correct behavior at x — > 0, when it goes to the Coulomb energy 
of an isolated nucleus, and at x — *■ 1, when it gives zero (electron and proton charges fully 
compensate each other). One can see that the Coulomb energy of the nucleus is screened by 
the electrons. This screening effect is very significant even at moderate electron densities. 
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For instance, at x = 10~ 3 (kp ~ 25 MeV) the screening effect is about 15%. As we shall see 
below, this change should lead to important modifications of the decay channels involving 
charged particles, such as a-decay and fission. 

The influence of the spherical electron cloud on the self-consistent proton potential is 
demonstrated in Fig. [3] for 208 Pb. We note that this nucleus, in the presence of electrons, 
is not stable anymore and would undergo electron capture. It is displayed here to demon- 
strate the effect of the attractive interaction between electrons and protons. The electron 
background leads to a downward shift of the potential but its structure remains basically 
unchanged within the nuclear volume. However, the principle difference from the pure 
Coulombic potential is evident at larger distances where the electric potential reaches zero 
(and zero derivative) at the boundary of the Wigner-Seiz cell. This is clearly seen in Fig. 
[3] for kp = 0.5 fm -1 and the cell radius Rc = 16.7 fm. This behavior is a result of the 
charge-neutrality condition which makes the electric field vanish outside the cell. We can 
conclude that at fixed Z and N the nuclear structure, in particular the single-particle level 
spacings, remains to good approximation unaltered by the presence of the electrons. As 
an example, we note that the proton rms radius in the heavy nucleus 240 Pu decreases by 
approximately 0.5% when adding electrons with kp = 0.5/m™ 1 . 

The Fermi energy of the protons also experiences a shift downward equal to the amount 
given by the additional electric potential. The neutron single-particle properties are only 
minimally altered, as is expected. The proton single-particle levels of 208 Pb are shown in Fig. 
H]for kp = 0.0 fm" 1 and kp = 0.5 fm" 1 . The downward shift of all levels is approximately 
10 MeV. We note that while the level density as such remains the same, the positioning of 
the levels with respect to the continuum threshold has changed. For kp = 0.5 fm , the 
number of proton states which correspond to bound orbitals has increased, and the density 
of states close to the continuum has increased. Furthermore, the proton separation energies 
have become larger. These changes in the single-particle spectra could affect, for example, 
nuclear reactions, in particular electron and neutron capture rates. 

The total charge density for the doubly-magic lead nucleus embedded in an electron gas 
in displayed in Fig. [51 The negative contribution due to the electrons and the downward 
shift in the interior of the nucleus are clearly visible. 
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FIG. 3: The proton (left) and neutron (right) single-particle potentials in 208 Pb for various electron 
Fermi momenta. The proton and neutron Fermi energies are shown as horizontal bars 

2. The deformed cell 

For spherical nuclei the spherical shape of the WS cell is the obvious choice. In this case, 
the electric field and its derivative vanish on the cell boundary so that different cells do not 
experience Coulomb interactions. However, when we investigate deformed nuclei, there are 
several possibilities of dealing with such a situation. One possibility is to employ again a 
spherical WS cell, with the deformed nucleus sitting in its center. This again corresponds 
to a constant electron background that is not affected very much by the presence of the 
nucleus. However, due to the deformed proton distribution, the quadrupole moment of the 
whole cell is nonzero in this case. This means that the neighboring cells will experience 
quadrupole-quadrupole interactions. We think that for the description of deformed nuclei 
it is more reasonable to use deformed cells too. Therefore, we consider axially deformed 
spheroidal cells with the excentricity determined by the condition of vanishing quadrupole 
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FIG. 4: (COLOR ONLINE): The proton single particle levels for k F = 0.0 fm~ l (left) k F = 
0.5 fm" 1 (right) in 208 Pb. The proton Fermi energies are shown as horizontal dashed lines 

moment. 

The quadrupole moment of a given charge distribution is defined as 

Q20 = \\[^ J d 3 xp(r)(2z 2 - r 2 ) (2.13) 

In our case p(f) = p p {f) — p e (f)- Thus, we adjust the shape of the electron spheroid 
such that Q20 = 0. For the electron Fermi momenta considered in this work the electron 
charge density is considerably smaller than the proton charge density at the center of the 
nucleus, and therefore the electron spheroid has larger axes. At the same excentricity, this 
would result in a larger negative quadrupole moment of the electron cloud as compared with 
the positive quadrupole moment of protons. Thus, to balance these two contributions, the 
electron spheroid must have a smaller excentricity than the proton one. 
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FIG. 5: The total charge density in a cell containing 208 Pb for various electron Fermi momenta 

The realization of the deformed WS cell is illustrated in Fig. [6] for the ground state and 
for a largely-deformed state of 240 Pu. The latter shape is chosen arbitrarily to illustrate the 
point concerning the shape of the WS cells. At the electron Fermi momentum kp = 0.5 frrT 1 
the equivalent spherical WS cell radius is 17.5 fm. In the ground state of Plutonium with 
a deformation of = 0.28, the electron cloud is almost spherical. At the deformation 
of fa = 1-9, the electron background clearly has a spheroidal shape. Still, its excentricity 
is much smaller than the one of the plutonium nucleus which is much more elongated. 
This shape corresponds to the nucleus behind the second barrier (see below). Note that 
this nucleus will eventually fission through the asymmetric barrier which is energetically 
favorable. However, the spheroidal shape used in our study can be used for axially-symmetric 
nuclei and, acoordingly, only for symmetric fission. Generalizations of this parametrizations 
would involve more complicated shapes allowing also for hexadecupole and octupole degrees 
of freedom. 

For the superheavy nucleus 292 120, the corresponding charge densities are shown in Fig 
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FIG. 6: (COLOR ONLINE): Total charge density of 240 Pu in the ground state (top) and at a 
deformation of 02 = 1.9 (bottom) for electron Fermi momentum kp = 0.5/m _1 



[3 Again, the (spherical) ground-state and a state at large deformation are displayed. This 



nucleus exhibits a semi-bubble 
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30] shape corresponding to a depletion of the baryon in 



its center that changes with deformation. 
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FIG. 7: (COLOR ONLINE): Total charge density of the superheavy nucleus 292 120 in the ground 
state (top) and at a deformation of = 1-5 (bottom) for electron Fermi momentum Hf = 0.5/m _1 
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FIG. 8: (COLOR ONLINE): The line of ^-stability (line in the middle, black), the proton drip-line 
(left-most line, blue) and the two-neutron drip-line (right-most line, red) for various electron fermi 
momenta. Magic numbers are indicated by grey lines. 

III. NUCLEAR CHART 



A. /3-stability line 

In this section we analyze how the nuclear /3-stability valley and the proton and neutron 
drip lines change with increasing electron Fermi momentum. As shown before, the proton 
single-particle potential experiences a constant downward shift. This downward shift of the 
proton potential and hence the gap between proton and neutron chemical potentials is the 
reason for the new /^-stability line. In the following, we investigate the possibility of certain 
decay modes of nuclei in the electron background. In particular, electron capture by nuclei 
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leads to a shift of the stability line to the neutron-rich side. 

In the presence of electrons the equilibrium condition with respect to weak decays (n — ► 
p + e~ + u e , p — > n + e + + u e ) reads 

Mn = A*p + /^e, (3.1) 

where \i n and fi p are neutron and proton chemical potentials calculated within the framework 
of the RMF model, i.e. 

fi N = J (m N + g a cr) 2 + p 2 FN + g^uj + g p p°r 3 , N = n,p (3.2) 

The electron chemical potential is simply given by 

fi e = yfk 2 F + m? e (3.3) 

where the contribution of the electrostatic potential has been neglected since it cancels out 
with the corresponding contribution from the proton in Eq. (13.11) . 

If /j, n > p, p + /i e holds, neutrons can decay into available proton states with the electron 
going on the top of the Fermi distribution. On the other hand, if \i n < fi p + /x e , protons 
can capture electrons from the Fermi distribution and occupy neutron states. Only when 
the condition in Eq. ( 13. II) is fulfilled, nuclear systems are stable with respect to these decay 
modes. Note that this relation becomes the well-known condition for /3-stability in vacuum, 
fi n fip, for h F = 0. 

For zero temperature, the proton and neutron chemical potentials coincide with their 
Fermi energies. We note that this correspondence is not uniquely defined anymore for 
systems with pairing, where the effective Fermi energy is determined to yield (N) = A, but 
can lie somewhere between the last level contributing with nonzero occupancy and the next 
one. Thus, in order to obtain the (N,Z) dependence of the stability line, instead of Eq. ( 13. lj) 
we employ the criterion 

where A = 1 MeV. 

In Fig. [HI the /3-stability lines are plotted for various values of kp. In the same plot 
we show also the proton and neutron drip lines (see below). As kp increases, the line of 
/3-stability is shifted more and more to the neutron-rich side, i.e. the presence of electrons 
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stabilizes neutron-rich nuclei. This shift is so strong already for hp = 0.1 fm~ l that /in- 
stability is reached only in the region of the two-neutron drip line. This means that the 
/3-equilibrium is reached only when free neutrons appear in the system. 

B. Proton- and neutron drip lines 

The neutron and proton drip lines are defined by the conditions \i v = m p and \i n = m n , 
respectively. It is clear that these nuclei are unstable with respect to weak decays so that 
the condition of /3-equilibrium, Eq. (13.11) . does not hold. In finite systems the dip lines 
can be defined in terms of separartion energies S p = B(N,Z) — B(N,Z — 1) and S n = 
B(N, Z) —B(N — 1, Z). Since in this work we are only calculating even-even system, we will 
consider instead two-proton and two-neutron drip lines. The two-proton drip-line is defined 
as the position where the two-proton separation energy, defined as 

S 2p (N, Z) = B(N, Z) - B(N, Z-2), (3.5) 

goes from positive to negative values. As an approximate relation, the single-proton sep- 
aration energies S P (N, Z) are related to the two-proton separation energies by S2 P (N, Z) « 
2 • S P (N, Z). The neutron drip line is the same for all cases shown in Fig. [8] since its position 
is not affected by the presence of the electrons. This follows from the fact that the neutron 
single-particle potential is not affected by the changes in the proton potential. Thus the 
neutron single-particle states are (almost completely) identical in all cases and the neutron 
drip line does not change. However, the results on the neutron drip line need to be taken 
with a grain of salt. Firstly, the calculations here are performed in spherical symmetry and 
deformation effects will wash out some of the strong shell effects appearing in these calcu- 
lations. Secondly, we employ BCS pairing which overestimates the coupling to continuum 
states and the pairing contribution of the loosely bound orbitals. For nuclei close to or at the 
neutron drip line, full Hartree-Fock-Bogolyubov calculations should be performed. However, 
since here we are not interested in the precise details of the position of the neutron drip line, 
the BCS calculations can give us useful information on the relation between the line of j3 
stability and the neutron drip region. 

The evolution of the proton drip-line as a function of the Fermi momentum is - in 
addition to Fig. |S]- displayed separately in Fig. As expected, this drip line shifts to more 
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FIG. 9: (COLOR ONLINE): The proton drip-line for various electron Fermi momenta. Magic 
numbers are indicated by grey lines. 

proton-rich nuclei as the electron Fermi momentum increases due to the additional binding 
of protons produced by the electron background. It is interesting that for kp = 0.5/m _1 , 
the proton drip line is very close to the N = Z condition. We can summarize our findings 
by stating that the region of nuclei between the conditions of proton drip and (3 stability 
has increased. 



IV. DECAY MODES 
A. a-decay 

In this section we investigate the influence of the electron background on a-decay of 
nuclei. The Q a value of the reaction is defined as 

Q a (N, Z) = B(N, Z) - B(N - 2, Z - 2) - B(2, 2) (4.1) 

and corresponds approximately, neglecting nuclear recoil, to the kinetic energy of the a 
particle leaving the nucleus. In the considered environment, the binding energies of mother 
and daughter nuclei as well as the binding energy of the a-particle are increased due to 
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FIG. 10: (COLOR ONLINE): Q a -values for k F = fin" 1 (top) and k F = 0.25 fin" 1 (bottom). 
The lower boundary corresponds to Q a = 0. Nuclei for which no convergent solution has been 
achieved have been left out. 

the attractive electromagnetic interaction of protons with the electron background. Fig. [10] 
demonstrates this effect. One can see that the overall trend is that the Q a values decrease. 
Also, the boundary for which Q a = shifts to more proton-rich nuclei. For Q a < 0, a-decay 
is not possible anymore and these systems are completely stabilized with respect to this 
decay mode. 

The calculation of the a-decay lifetimes in the presence of electrons is not completely 
trivial. There are two competing effects. As seen above, the Q a value is lowered due to 
the increased nuclear binding caused by electrons, which alone (for the same barrier) could 
increase the lifetime. We have also seen that the presence of the electrons leads to screening 
and modification of the Coulomb potential, see Eq. (12.121) . Therefore, we need also to take 
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FIG. 11: a-potential for a-decay of the nucleus 218 U. 

into account the change of the barrier through which the a particle has to penetrate. To 
quantify this effect we use a simple one-parameter model of Ref. [3JJ which expresses the 
a-particle potential V a as 

V a = 2V p + 2V n , (4.2) 

where V p and V n are, respectively, proton and neutron single-particle potentials. The a 
particle is treated as a boson. The half-life is written as T1/2 = In2/A, where the decay 
constant A is parametrized as A = c P. Here the pre-formation factor and the knocking 
frequency in the Gamov picture are absorbed into one parameter c, which is adjusted to 
known data l . The probability for transmission through the barrier, P = e s , is calculated 
within the WKB approximation 

S = -2 



R2 dr j MV a {r) -Q c 



(4.3) 



where R\ and R2 are the tur ning points of the barrier. 

As demonstrated in Ref. 3J|, this model performs better than the 4-parameter Viola- 
Seaborg systematics and additionally incorporates isovector trends, i.e. the a lifetime for 



1 A model variant that computes the knocking frequency numerically 31 1 yields quite similar results 
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FIG. 12: The evolution of Q a values (left) and a decay half-lifes (right) for the nucleus 218 U as 
a function of the electron Fermi momentum. In the right panel, the full squares correspond to 
calculations within the model of Ref. 3l|, and open squares give results for life times multiplied 



with the prefactor of Eq. (|4.6p . The lines are drawn to guide the eye 

a nucleus with given charge number possesses a dependence on the neutron number due to 
the self-consistent interrelation between proton and neutron single-particle potentials. 

We employ this model here also for finite electron density where the potentials and the 
Q a values are taken from the RMF calculations in the electron background. As an example 
we present results for the nucleus 218 U with the magic neutron number N = 126. We have 
checked in axially-symmetric calculations that both this nucleus as well as its daughter are 
spherical, thus no deformation effects enter the results on a-decay. 

As seen in Fig. [TTJ the a-particle potential is lowered as a function of increasing electron 
Fermi momentum, making it easier for the a particle to escape the nucleus. Together with 
the Q a value this leads to an overall decreasing of a half lives, see Fig [12j This holds true 
at low electron densities. However, this trend changes in the vicinity of the point where 
the Q a values turn negative at some electron Fermi momentum kp. At kp > kf lt , a-decay 
becomes forbidden and nuclei are stable with respect to this decay mode. In the Uranium 
isotope discussed here, kp lt = 0.24 jmT 1 . However, the potential barrier experienced by 
the a particle at k c p lt = 0.24 fmT 1 is still « 13 MeV high. If the barrier would result from 
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a pure nuclear charge, the lifetime for Q a ps would increase by a huge amount compared 
to situations with a finite Q a value. The reason is that in vacuum the Coulomb barrier 
goes to zero asymptotically as 1/r, thus at Q a ~ the integration distance for the WKB 
formula is very large. By using a Taylor expansion of the Gamov factor for a 1/r potential 

— 1/2 

one can find the approximate relation lnri/2 ~ Q a , which leads to an infinite lifetime 
at Q a — > 0. In the case of electron background, however, the Coulomb potential is fully 
screened for distances larger than the Wigner-Seitz radius, which decreases with increasing 
electron Fermi momentum. Thus, this approximate relation is not applicable anymore. We 
expect that the WKB-based calculation underestimates the upwards trend shown in Fig. [121 



(right), since close to the threshold this method is not accurate 32]. 

This situation can be illustrated by a text-book 1-D tunneling problem. A particle coming 
from the negative x direction_encounters a rectangular barrier of height V between x = 
and x = d. It can be shown 



33] that the transmission coefficient T is given by 



T 



1 + V - S inh 2 f ^-^ l " (4 4) 



where E and m are the energy and the mass of the tunneling particle, respectively. 
For typical cases of E ~ V/2 and ^j2m{V — Wjd/h >> 1 (corresponding to a broad 
barrier compared to the de Broglie wave length of the tunneling particle), we have 
sinh( v /2m(F - E)d/H) w (1/4) exp(2^/2m(V - E)d/h) >> 1. Since the pre-exponential 
is of order 1, we can neglect the first term in brackets of Eq. (14.41) and obtain 

T w exp{-2/h v / 2m(V - E)d] (4.5) 

This is exactly what follows from the WKB approximation of Eq. (14.31) . For electron Fermi 
momenta close to kp lt , however, we have E — ► 0, and hence the prefactor diverges, i.e., 

V 2 V 

oo, E -> (4.6) 



V 2 - (2E - V) 2 AE 

This gives additional suppression of the tunneling probability as compared to the one ob- 
tained in the WKB approximation. Thus, the trend to longer lifetimes can be expected 
to set in for even smaller values of kp compared to our calculations using the WKB ap- 
proximation (Fig. [121 full squares). This means that there is a k^ m for which the lifetime 
reaches its minimum, and for kp > /c™ m the lifetimes increases until the nucleus becomes 
stable with respect to this decay mode for kp > kp lt . To illustrate the trend we have made 
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estimates of the lifetimes obtained by multiplying Eq. (14. 5 j) with the prefactor of Eq. (14. 6 p 
for E = Q a . Results are shown in the right part of Fig. [12] with open squares. In this 
case the reversal of the downward trend in half lives occurs earlier and more dramatically 
compared to the previous calculations. For the electron Fermi momentum kp = 0.1 fm _1 
the prefactor has the value 1.53, while for kp = 0.2405 fm -1 it reaches the value 1255. A 
further modification would arise from the fact that the exponential function in Eq. (14. 4p 
would have to be replaced by the sinh function of Eq. (14. 5 p which could yield even more 
stabilization with respect to a decay. 

The linear decrease of the Q a values shown in Fig. [12] can be understood if we consider 
the electric potential caused by a spherical electron cloud given in Eq. 12.101 The potential 
energy V pot of a point charge located in the middle of the cloud is inversely proportional to 
the WS cell radius Rc- Taking into account that Rc oc 1/kp, Eq. 12.61 we obtain V pot oc kp, 
i.e., a point-like nucleus gains potential energy proportional to kp. If Z denotes the charge 
number of the mother nucleus and Z\ and Z 2 denote the charge numbers of the daugther 
products [Z\ = 2 and Zi = Z — 2 for a decay), The change in the Q value stemming from 
the additional binding due to the electron cloud is given by 

SQ = ^ (^) 1/3 e*k F [Z^ - ZT ~ ZT] (4-7) 

Thus the Q a value which depends on the energies of the mother and daughter nuclei, de- 
creases linearly with kp. We also see that since Q oc Z 5 ^ 3 the largest contribution is coming 
from the mother nucleus. 

B. Spontaneous fission 

1. Deformed calculations 

In order to investigate the influence of the electron background on the fission barrier of 
heavy and superheavy nuclei, we compute the energy of the system as a function of deforma- 
tion. The cuts through the potential energy surface (PES) are calculated in axial symmetry 
for reflection-symmetric shapes using a constraint on the total quadrupole moment Q20 of 
the nucleus. This is achieved by adding — AQ20 to the Hamilton operator and minimizing 
(H — XQ2o)- All other multipole moments that are allowed by the symmetry of the calcu- 
lation are not constrained and adjust themselves to the solution of minimal energy. So, in 
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contrast to the macroscopic-microscopic approach 34(, we are not operating in a limited 
deformation space. On the other hand, we are not necessarily exactly following the gradient 
in the multidimensional PES which would only be the case if the fission valley would be 



parallel to the Q20 direction 



35|. 



At each quadrupole deformation along the fission barrier, the electron clouds excentricity 
and spatial extension is iteratively adjusted to match two constraints: 1) the total charge 
of the deformed WS cell has to be zero, and 2) the quadrupole moment of its charge dis- 
tribution (including contributions from both protons and electrons) should vanish. These 
two constraints can be realized by varying two parameters characterizing the electron cloud, 
for example the lengths of long half-axis and the excentricity. As discussed above, these 
conditions minimize the mutual interactions of WS cells. 



2. Fission barriers 

In fission studies an important role is played by the energy of the system as a function 
of the deformation parameter Our analysis shows that the shape of the electron cloud 
has some effect on the energy. While the spherical cloud does not produce any changes, the 
deformed cloud leads to visible effects. 

Important key quantities related to fission are the width and the height of the fission 
barrier. The influence of the electrons on the fission barrier is a macroscopic effect resulting 
from interaction of the protons with the Coulomb potential produced by the electrons. The 
electrostatic repulsion between the protons is weakened by the presence of the negative 
charge background. This change leads to the increase of barrier height and isomer energy. 
Physically speaking, the electron background tends to slightly stabilize the system with 
respect to symmetric spontaneous fission and to increase the excitation energy of the isomeric 
state. We expect that asymmetric fission will be altered in a similar way, although our 
consideration does not allow octupole shapes. We note that the potential energy surface is 
determined by the dependence of the energy on the deformation. Especially for superheavy 
systems, the liquid drop barrier vanishes, and their stabilization results from shell effects. 

The trends discussed above are clearly seen in Fig. [13] which displays the fission barrier of 
Plutonium for different choices of the electron density in the spheroidal configuration. The 
inner barrier increases with increasing hp leading to a stabilization effect towards sponta- 
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FIG. 13: Fission barrier of 240 Pu for various electron densities in spheroidal configurations. 



neous fission. Furtheron, the properties of the shape isomer (second minimum) are altered 
due to the electrons background too. For kp = 0.5 fm we see an increase of the inner fission 
barrier by approx. 0.3 MeV and the energy of the isomeric state increases by approximately 
1 MeV. Note that the ground-state binding energy of 240 Pu with NL3 is 1814 MeV, the pres- 
ence of the electrons increases it to the value of 2248 MeV (the latter does not include the 
kinetic energy of the electrons). As another example, we show the influence of the electron 
background on the fission barrier of 256 Fm in Fig. [T4] The presence of the electrons leads 
to a slight increase of the first barrier. The effect is more dramatic for the isomeric state, 
where an increase in energy of 2 MeV is found at kp = 0.5/m -1 . Note that the low energy 
of the second minimum, which is below the ground state, is a (unrealistic) prediction which 



36|]. 



has also been found in a systematic study for heavy and superheavy nuclei 

Fig. [15] demonstrates the influence of the electron background on the superheavy nucleus 
292 120. In this case, the electron background leads to a slight decrease of the width and 
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FIG. 14: Fission barrier of 256 Fm in for various electron densities in spheroidal configurations. 

height of the inner barrier calculated with the spheroidal electron cloud. But the height of 
the second minimum has increased by about 0.7 MeV. Note, however, that the inclusion 
of reflection-asymmetric degrees of freedom in these superheavy systems make the second 
barrier vanish [36( . It is interesting that the potential energy surface is affected in a different 
way compared to the case of Plutonium. The reasons are related to the fact that the barrier 
starts out from a spherical minimum, and that the charge number in this nucleus is much 
larger so that the total contribution of the Coulomb energy to the total energy is strongly 
increased compared to plutonium. 



We can speculate about even heavier systems. In Ref. [37(|, the shell corrections for 
super- and hyperheavy nuclei were calculated within the framework of self-consistent mean- 
field models. It was found that the familiar concept of strongly-pronounced magic numbers 
dissolves due to the high level- density, and rather a broad region of shell-stabilization can 
be found. We choose the hyperheavy nucleus 462 1 54, which is located in such a region of 
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FIG. 15: Fission barrier of the superheavy nucleus 120 for various electron densities in spheroidal 
configurations. 

shell-stabilization (for the RMF force NL3), as an example. Its fission barrier can be seen 
in Fig. [T6j The strong shell corrections lead to a spherical minimum and to a barrier which 
is which is about 6.5 MeV high, and almost 1 MeV larger in the presence of electrons. 

A word of caution is in order. For both this hyperheavy nucleus and the superheavy 
system 292 1 20 discussed before, the calculation of the axial barrier gives us only partial 
information on the principal impact of the electron background on the fission barrier. The 
actual fission path, however, might deviate quite strongly from this path, and for these heavy 
systems fissioning through the triaxial plane is rather probable. 
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FIG. 16: Fission barrier of the superheavy nucleus 462 154 without electrons and for \zf = 0.5 fm _1 
3. Fission Q values 

As discussed in Section III B 1[ due to the screening of the electrons the Coulomb energy 
of nuclei is lowered as compared with the vacuum case. This will lead to the reduction of 
the Q value for spontaneous fission, in full analogy to the case of a decay. Fig. [17] shows 
the Q value for symmetric fission of 256 Fm as a function of the electron Fermi momentum. 
One can clearly see this trend and we think it will be present for asymmetric fission too. To 
evaluate the fission probability one should perform detailed calculations of the fission barrier 
for the case of a two-center nuclear shape. This work is in progress now. At this stage we 
can only point out that due to the screening of the Coulomb potential the fission barrier is 
altered and the Q value will decrease. We expect that the latter trend is more important 
so that the net effect is suppression of spontaneous fission. This follows also from simple 
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FIG. 17: Q value for symmetric fission of 256 Fm as a function of the electron Fermi momentum 
estimates based on the liquid-drop model. Indeed, let us consider the fissility parameter 



Z 2 2a 



(4.8) 



A a c 

where as ~ 18 MeV and ac ~ 0.72 MeV are, respectively the surface and Coulomb 
coefficients in the Weizsacker formula. If we take into account the electron screening, the 
Coulomb energy is reduced by the factor c, Eq. (I2.12p . Therefore the modified fissility 
parameter is ^ = — . Since c(x) < 1, this means that in the presence of electrons the region 
of spontaneous fission moves to heavier nuclei. For instance, for kp = 0.25/m -1 , the fissility 
parameter is 80 instead of 50 in vacuum. 



V. CONCLUSIONS AND OUTLOOK 



We have studied the properties of atomic nuclei embedded in an electron gas as occuring 
in, e.g., neutron star crusts and supernova explosions. Nuclear structure calculations have 
been performed within the relativistic mean-field approach employing the force NL3. The 
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calculations have been performed within the Wigner Seitz cell approximation in coordinate- 
space. A fermi-type distribution with a smooth surface has been employed in order to get 
rid of artifacts stemming from the discrete grid spacing. For spherical systems, a spherical 
Wigner-Seitz cell has been used. We have discussed an implementation of spheroidal Wigner- 
Seitz cells for axially-deformed nuclei. We employ the criterion that in addition to the total 
charge also the total electric quadrupole moment of a cell should vanish. This leads to an 
iterarively self-adjusting Wigner-Seitz cell for each given deformation of the nucleus. 

The presence of electrons leads to effects which have rather macroscopic character. We did 
not find any significant changes in the single-particle levels regarding their relative positions. 
Since the electron background is constant over the nuclear volume - by construction - it 
leads predominantly to a constant downward shift of the proton potential. 

The electron gas alters the stability condition for (3 decay and shifts the stability line 
to more neutron-rich nuclei. Electron capture becomes a relevant process, the transforma- 
tion of protons into neutrons favors large isospin in nuclei. Furthermore, the two-proton 
dripline is shifted to more proton-rich nuclei since the protons gain additional binding due 
to the attractive interaction with the electrons. The neutron drip line is not altered by the 
presence of the electrons since the neutron single-particle potential remains (to a very good 
approximation) unaffected. 

We have found that in a dense electron background decay modes such as a-decay and 
spontaneous fission are suppressed. Our calculations show that the a decay half lives de- 
crease a a function of the electron Fermi momentum fcp until they get stabilized for larger 
k F . A general trend is that the Q values of those decay modes decrease with increasing 
electron density. 

We have also calculated the potential energy surface as a function of the electron Fermi 
momentum within the framework of the self-adjusting axial Wigner Seitz cell. We have 
found that only for very large electron Fermi momenta a change is occuring for the inner 
axial barrier and the energetic position of the isomeric state. This could hint again at 
stabilization of the fission mode for high electron densities. 

Overall, the electron gas broadens the part of the nuclear chart lying between the proton 
drip line and the valley of (3- stability. Stabilizing effects with respect to a decay and 
spontaneous fission occur for large electron Fermi momenta. Hence, in extreme astrophysical 
environments, the production of very exotic and superheavy nuclei could become possible. 
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This might happen during the r-process (rapid neutron capture) and the rp-process (rapid 
proton capture) when the electron background prevents the heavy (and superheavy) nuclei 
from fast decay by spontaneous fission or alpha decay. As a consequence, the nuclei, which 
otherwise would be unstable, can provide a bridge to the island of superheavy elements. We 
are planning to study this possibility in the future. If long-lived or even stable superheavy 
nuclei exist, they could be created in such an environment and later ejected into space. Then 
one can try to search for such superheavy nuclei in cosmic rays. Another natural step in our 
investigations is to add free neutrons and study nuclei in the environment of free electrons 
and neutrons. Work in this direction is in progress. 
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Appendix 

A. Role of Coulomb exchange 

In relativistic models such as used in this study, the Coulomb interaction is usually 
calculated within the Hartree approximation, i.e. 



£Toui omb = e 2 -J J tfVdV ' (5.1) 

where p p (f) is the proton density. This relates to the fact that the RMF model traditionally 
is formulated as a Hartree theory, and no explicit exchange terms in the effective nucleon- 
nucleon interaction are taken into account. In contrast to effective field theories with contact 
interactions, exchange terms involving finite-range meson fields are computationally expen- 
sive, and so far no significant improvement over Hartree-like formulations has been achieved. 



The role of 4-fermion interactions has been studied in Ref. [38| for the point-coupling variant 
of the RMF model, the RMF-PC model. It has been found that, while the formal struc- 
ture of the model remains the same and only a redefinition of 4-fermion coupling constants 
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occurs, the interpretation of the various terms becomes quite different. For higher-order 



39] yield a 



point-coupling terms and derivative terms, however, the Fierz transformations 
large number of terms which do not improve the model 401 ]. 

In the Skyrme-Hartree-Fock approach, the Coulomb exchange term is usally included in 
Slater approximation 4=11 ] : 

Coulomb = -/ (I) V3 / d 3 r [p p (r)f 3 (5.2) 

In our calculations, we replace theproton density by the absolute value of the charge density, 
i.e., p p (f) by \p ch {r)\ = \p P (r) + p e (r)\- 

For nuclei close to stability, it can be shown that the Coulomb exchange effect can be 
absorbed by a refit of the mean- field parameters 42| . It remains to be seen, however, if this 
still holds true for exotic nuclei and nuclei close to the proton drip-line, as well as nuclear 
systems studied here. As demonstrated in Fig. [18], the Coulomb exchange interaction both 
in vacuum and at high electron chemical potential yields a rather small contribution to the 
potential, and Fermi energies for the two cases differ by 0.4 MeV only. Thus the properties of 
nuclei in an electron gas are not sensitive to the exclusion or inclusion of (Slater-) Coulomb 
exchange in the calculation. 
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FIG. 18: The proton single-particle potential for 218 U in the absence of electrons (left) and for 
kp = 0.5/m _1 (right) involving only the direct Coulomb term (H) and involving both direct term 
and exchange term in Slater approximation (HF). Also shown are the proton fermi energies for 
both cases. 
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